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Abstract. M dwarfs are known to flare on timescales from minutes to hours, with 
flux increases of several magnitudes in the blue/near-UV. These frequent, powerful 
events, which are caused by magnetic reconnection, will have a strong observational 
signature in large, time-domain surveys. The radiation and particle fluxes from flares 
may also exert a significant influence on the atmospheres of orbiting planets, and affect 
their habitability. We present a statistical model of flaring M dwarfs in the Galaxy that 
allows us to predict the observed flare rate along a given line of sight for a particular sur- 
vey depth and cadence. The parameters that enter the model are the Galactic structure, 
the distribution of magnetically active and inactive M dwarfs, and the flare frequency 
distribution (FFD) of both populations. The FFD is a function of spectral type, activ- 
ity, and Galactic height. Although inactive M dwarfs make up the majority of stars in 
a magnitude-limited survey, the FFD of inactive stars is very poorly constrained. We 
have organized a flare monitoring campaign comprising hundreds of hours of new ob- 
servations from both the ground and space to better constrain flare rates. Incorporating 
the new observations into our model provides more accurate predictions of stellar vari- 
ability caused by flares on M dwarfs. We pay particular attention to the likelihood of 
flares appearing as optical transients (i.e., host star not seen in quiescent data). 



1. Introduction 

Stellar flares are powered by magnetic reconnection events. Alth ough most last for onl y 
a few minutes, they have been observed to last for up to 8 hours ( Kowalski et alj|2010h . 



Their strongest observational signal is a flux increase in the blue/near-UV, which can 
be several magnitudes for the biggest flai^es. Because M dwarfs account for such a large 
fraction of the stars in the Galaxy, photometric variability on these stars will ha ve an 
significant signature in time- domain survey s such as Pan-STARRS (Kaiser 2004), PTF 
(iLaw et al.ii2009) . and LSST (ILSST Science Collaborations: Paul A. Abell 2009) . 



M dwarfs are an excellent place to look for extrasolar planets because planets have 
relatively large RV effects and transit depths on low-mass stars. Many planets, including 



relatively large KV effects and transit depths on low-mass stars. Many planets, including 
several super-Earths, have been discovered o rbiting M dwarfs (e.g., lUdry et aLll2007l : 



ICharbonneau et al]|2009l : ICorreia et alj|2010l) . However, magnetic activity on M dwarfs 



can hamper planet search es, since starspots and flares introduce both photometric nois e 
and radial velocity jitter (IWrightll2005| - iBasri et ai]|2010l : lL6pez-Santiago et alj|201(]h . 



At the same time, the intense flux of high energy photons and particles caused by flares 
may have a significant impac t on the atmospheres of Earth-like plan ets in the close-in 
habitable zones of M dwarfs dBuccino et allboOTl : ISegura et alJlioiOh . 
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Our group has used the serendipitous obse rvations of flares in surveys to make 
estimates of the flare rate. Kowalski et al. (I2009h searched for flare s in th e ~80 photo 



metric epochs from the SDSS Stripe 82 survey (llvezic et al.ll2007n . while iHilton et al 



CpiO) identified flares in the sample of M dwarfs with SDSS spectroscopy (IWest et al 
120081) . These studies confirmed that while the flare rate increases with later spectral 
subtype, the flare energy on earlier type stars was higher. We also found that the flare 
stars were preferentially closer to the Galactic midplane, which we interpret as an age 
effect, with younger stars flaring more frequently. 

Although the occurrence of an individual flare cannot be predicted, dedicated pho- 
tometric monitoring campaigns of individual stars can yield measurements of the fre- 
quency of flares of various energies, also called the flare frequency distribution (FFD). 



Mofi^ettl(fT974.') and iLacv et all (fl97^ obtained FFDs on several of the most well-known 



and most prolific flare stars in the Solar neighborhood, nearly all of which are active 
mid-type M dwarfs. We have extended these measurements of the FFD to inactive stars, 
as well as both early- and late-type M dwarfs. 



2. Observations 

Measuring the FFD of a star requires time-resolved photometric monitoring over the 
course of many flares. We have used four telescopes (described in Table 1) to obtain 
almost 500 hours of flare monitoring of our targets. 



Table 1 . The telescopes used to collect flare-monitoring data 



Location 


Size 


Operator 


Filters 


Apache Point Observatory 


3.5m 


Astrophysical Research Consortium 


B 


Apache Point Observatory 


1.0m 


New Mexico State University 


U 


Manastash Ridge Observatory 


0.8m 


University of Washington 


u,g 


Apache Point Observatory 


0.5m 


Astrophysical Research Consortium 


u 



We chose bright targets in order to have higher time resolution while maintaining 
high signal-to-noise. The cadence of our observations was typically about ten seconds, 
except for the NSMU 1.0m telescope, which has a longer readout time and a cadence of 
approximately 60 seconds. The data reduction and aperture photometry was done with 
standard IRAF tasks. Because we are only concerned about variation in the lightcurve, 
we used differential photometry. We use the equivalent width of Hor to classify each of 
our tar gets as either magnetically active or inactive. Following the criteria of West et all 
(l2008h . stars with Ha equivalent width of more than lA in emission are considered 
active. 



3. Flare Frequency Distributions 

In order to derive the flare frequency distributions for each bin of spectral type and 
activity level, we first need to identify flares, calculate the flare energies, transform 
each of these flare energies to a common filter, and then combine the observations from 
different nights and different stars. 
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3.1. Flare Identification and Energies 

For each night of observations, we first define a period of quiescence. This is identified 
by eye and is at least an hour in duration. We use this quiescent period to define the 
quiescent flux (taken as the median flux value) and the precision of our data, taken as 
the standard deviation of the data points during the quiescent period. We define a flare 
as having occurred if there are at least three consecutive measurements more than 2.5 cr 
above the quiescent value, with at least one of those measurements being > 5 cr above 
the quiescent value. The flare start and stop times are identified as when a running 
average of several measurements is equal to the quiescent value. Since this procedure 
requires the flux to return to the quiescent value before the flare is determined to have 
ended, groups of flares or sub peaks within a flar e are typically id e ntified as one flare . 

We follow the method of Gershberg (19721); iMoffetd (1 19741) : iLacv et all (Il976h to 
find the flare energies in the observed filter. We first calculate the flare flux as 



I fit) = 1 



(1) 



where Iq is the quiescent flux. The equivalent duration, P, is simply the time integral 
between the time of flare start and stop (described above). 



P = J hit) 

The total flare energy in filter x is 



dt 



(2) 



(3) 



where qx is the quiescent flare luminosity in filter x. The quiescent flare luminosity is 
easily calculated from the known absolute magnitude of the star. The abs olute mag- 
nitud e s come from the Palomar/MSU Nea r by Star Spectros copic Survey (iReid et al.1 
[1996'; 'Hawley et al.|ll997|; iGizis et al.ll2002l : iReid et al.ll200 2'). or the Gleise catalogue 
(IGliese 1969) . We transform the energies of al l flares into SPSS u using our own si- 
multaneous observations and a relationship from lLacy et al.l ( T976l) . The total observing 
time, number of stars, and number of flares in each bin is given in Table 2. 



Table 2. Flare Monitoring Targets 

Bin Total (hrs) # of Flares # Stars/Bin 

Inactive MO-2 201.05 5 9 

Inactive M3-5 107.96 4 3 

Active M3-5 123.15 43 7 

Active M6-9 60.27 29 4 



Total 492.43 hours 81 flai-es 23 stai's 



3.2. Constructing the FFD 

The FFD is a cumulative distribution of flare energies such that the FFD evaluated at 
energy E is the number of flares per unit time with energy greater than E. To construct 
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the FFD, we need not only the energy of each flare, but also the duration of observations 
during which a flare of that energy is able to be detected. Because conditions change 
from night to night, and because we must bin together stars of different magnitudes 
with different comparison stars, the minimum detectable flaie energy is different for 
each night on each star. 

For each bin, the FFD is constructed from the hst of flare energies in that bin and 
the duration (t) and minimum observable flare energy (£''"'") for each observing period 
for stars in that bin. The differential flare frequency is the inverse of the sum of all 
observing periods whose minimum observable energy is less than the flare energy. 

^ ma x 

FF{Ei) = 11 Yj^i (4) 

i 

where £"""" is the minimum observable flare energy on a night with duration r,, and 
^'"max is the night with the largest threshold for detecting a flare. The cumulative FFD 
is simply the cumulative sum of the flare frequencies. 

The cumulative FFDs for each o f our bins are shown in Figure 1. Previous studies 
(e.g., iLacy et al]|l976l : lAudard et"aD ,2000) have described the FFD as a linear relation 



between the log of the flare energy and the log of the frequency 

log V = a + y6 log ff/ (5) 

Overplotted in Figure 1 is our best fit for each FFD. 

The FFD for the active M3-M5 stars is consistent with the measured FFDs of 
several of t he mo st well-known flare stars (all active and between M3 and M5) from 



Lacy et al.l (119761) . The FFDs for the other three bins represent the first measurement for 
these bins. We note that inactive stars flare much less frequently than active stars, and 
with less energy. For the inactive stars, flares on earlier type stars are more energetic 
than on later type stars. This trend is also apparen t for the active s tars, if the FFD 



of the only measured active early M star, Au Mic ( Lacy et all 1 19761 not shown), is 
representative of all active early M stars. Although it may be that flares on earlier type 
stars are intrinsically more energetic, this is at least partially a selection effect, since a 
flare on an earlier type star will cause a smaller change in apparent magnitude than a 
flare of equal energy on a later type star, the so called "contrast effect." 



4. Galactic Flare Model 



We would like to predict the number and size of flares that will be seen in time domain 
surveys. To do this, we create a model of flaring in the Galaxy. The first step in model- 
ing the Galactic flare rate is r nodehng the light curves of individual sta r s. We adopt the 
classi c flare shape (see e.g.. lMoffettlll974l : Hawley & Petterseiilll99ll : [Kowalski et al. 



I2010h of an impulsive phase with a linear flux rise and partial linear decrease, followed 
by a gradual phase consisting of an exponential flux decay. We insert flares into a light 
curve by drawing from the measured FFD for a particular type of star. Using the model 
light curves, we calculate the fraction of time a star is seen at increased brightness. The 
black line in Figure 2 (top panel) shows the FFD for the active M3-M5 bin, while the 
black line in the bottom panel shows the fraction of time a hypothetical dM3.5e star 
is seen at increased brightness from flares. The uncertainty in the measurement of the 
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Figure 1 . Using nearly 500 hours of flare monitoring observations, we calculate 
the FFDs of four bins of differing spectral type and activity status are shown here, 
along with the best fit model. We find that inactive stars flare less frequently than 
active stars. 



FFD (shaded gray) is dominated by our choice of the flare definition, especially the 
start and stop times, and the assumption that all stars within a particular bin have the 
same intrinsic FFD. Dilferent input FFDs result in significant changes in the amount 
of time spent at increased brightness. In principle, we can use time domain surveys, 
which provide measurements of the fraction of time spent at increased brightness, to 
determine which FFD is the best fit to the survey data. 

We have created a model of flaring on M dwarfs in the Galaxy to predict the 
number and magnitude of flares that will be observed in time domain surveys, with 
a particular focus on L SST. We generat e a model Galaxy with thin disk, thick disk, 
and halo components ( Juric et al. |2010[) . The scale heights, lengt hs, and densities 



of these com ponents have been determined from SDSS star counts dJuric et al.l 12008 
Bochanski et al.. 2010) . The active fr actions of each sp ectral subtype as a function of 



Galactic height are extrapolated from lWest et all (12008 

For a given pointing and field of view size, the Galactic model determines the 
number, color, distance, and activity level of M dwarfs along that line of sight. The 
activity status and spectral type determine the FFD used to generate a light curve for 
each star in the model. 

To verify t he utility of this model, we compared our model results to the mea- 
sured flare rates (iKowalski et al.i r2009) from the SDSS Stripe 82 photometric catalogue 
(llvezic et aLll2007h . Each of the -50,000 M dwarfs in Stripe 82 was ob served -80 
times, with 271 epochs meeting stringent flare criteria dKowalski et al.ll2009h . The mag- 
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Figure 2. Top; The measured FFD of an active M3.5, with the uncertainty shown 
as shaded gray. Bottom: The fraction of time an active M3.5 spends at increased 
brightness. Stars with different FFDs spend vastly different amounts of time at highly 
enhanced magnitudes. Future studies may be able to use the large number of flare 
observations from time domain surveys to determine the FFD more precisely. 



nitude distribution of these flares is shown in Figure 3 (red line). Adopting the same 
number of epochs, field of view, flare criteria, and survey depth of the Stripe 82 data 
in our flare model, we predict 283 flares, which are overplotted in Figure 3 (blue line). 
Our model not only matches the total number of flares, but independently reproduces 
the distribution of flare magnitudes remarkably well. 
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Figure 3. The magn itude distribution of flares measured in SDSS Stripe 82 from 
iKowalski et alJ (l2009h (red line) and the model prediction using the same survey 
parameters (blue line). 



4.1. Predicting Flare Rates in Time Domain Surveys 

The good agreement between our flare model and the measured Stripe 82 flare rates 
gives us confidence that we can predict the frequency of flares in future time domain 
surveys. However, there are three areas where limited data may aff"ect the model pre- 
dictions. These limitations are especially important to the prediction of transients — 
flares that are bright enough to make a star briefly visible that is otherwise too faint to 
be seen in the survey. 

The first uncertainty is the active fraction of stars as a function of spectral type 



and distance from the Galactic midplane. The IWest et all (2008) results are based on 
SDSS spectroscopic data, which are limited to bright objects. The active fraction at late 
spectral types and at far distances, which will certainly be visible to LSST, are currently 
unconstrained, and will likely remain so in the near-term. The second reason is that a 
study of flares in the SDSS spectroscopic survey finds that flare sta rs are preferentially 



closer to the Galactic plane than the active stars ( Hilton et al.ll201Cll) . indicating that the 



FFD of active stars may be a function of distance from the Galactic plane. The third 
limitation to our model inputs is in the rate of flares at the high energy, low frequency 
end of the distribution. Because large flares happen infrequently, we have limited infor- 
mation on the frequency of the highest energy flares. We don't know if there is a break 
in the distribution, or simply a cut-off". 

We have used the expected survey parameters for LSST to predict the number of 
flares seen in a single 15-second w-band exposure. Figure 4 shows the number of flares 
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seen as a function of Galactic latitude. The decrease with increasing latitude is primarily 
caused by the smaller number of stars in sightlines looking out of the Galactic plane. 
As expected, there are more small flares than large flares. The error bars are calculated 
using the uncertainties in the FFDs (see Figure 2). The number of flares seen in a 



single night of LSST observing will be larger than the entirety of the iKowalski et al. 
(12009) SDSS Stripe 82 results. In the future, we can determine the FFD of stars to a 
much greater precision, and for much smaller spectral type and Galactic height bins, by 
varying our input FFDs in the model to fit the large number of LSST flare observations. 




40 60 
Galactic Latitude [degrees] 



Figure 4. The predicted number of flares that will be seen in a single LSST M-band 
exposure. The decrease with Galactic latitude is caused by a decrease in the number 
of M dwarfs in the sightline. 

Astronomers interested in exotic transients can neglect most stellar flares, since 
the source of the flare is known from the quiescent photometry. However, it is certainly 
possible for an M dwarf to have such a large flare that it will appear as an optical 
transient. In this case, the star is too faint to be visible in quiescence in any filter, but is 
briefly bright enough that it is seen during the flare, and will thus appear as a blue point 
source with no observed host. To meet this criteria, a flare must have a large change in 
magnitude: 

Am > {Ug - Zq) + {Zq - ZUm) + {ZUm " "/mi) = Uq - UUm (6) 

where Zq < Zum (7) 

where the subscript q represents quiescent measurements and the subscript Urn repre- 
sents the survey limiting magnitude in that filter. The first term is the u - z color of the 
star, which is several magnitudes for late M dwarfs. The second term means that stars 
close to the survey limit more likely to appear as transients, and the final term is the 
ditference between the survey limits in the two filters. Note that although we use u and 
z filters in this example, we can make similar predictions for any other filter. 

Using the LSST survey criteria, we predict that the probability of M dwarf flares 
appearing as optical transients in a single w-band exposure is ~10% at low Galactic 
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latitude, decreasing to just a few percent at high latitudes. Only the largest flares will 
meet the transient criteria, and as mentioned above, the frequency of very large flares is 
not well measured. The error bars are therefore quite large. 



5. Conclusions 

We have collected nearly 500 hours of flare monitoring observations to make the first 
measurements of the FFD of inactive early- and mid-type M dwarfs and for active late- 
type M dwarfs. The FFD is best represented as a linear relation between the log of 
the quiescent flare energy and the log of the frequency. We find that flares on earlier 
type stars are more energetic than on later type stars, and that inactive stars flare less 
frequently than active stars. 

We incorporate the measured FFDs into our Galactic M dwarf flare model, which 
can be used to predict the number and magnitude of flares that will be seen in any 
survey. To make this prediction, we also adopt a flare shape, an active fraction (which is 
a function of spectral type and distance from the Galactic plane), a luminosity function, 
and a Galactic model. We draw from the appropriate FFD to generate a light curve for 
each star along a given sightline. Using the survey parameters from SDSS Stripe 82, 
our model was able to reproduce bo th the number and the magnitude distribution of 



flares found bv lKowalski et al.l (l2009h . 



Applying our model to the survey parameters of LSST predicts tens of flares > 0.1 
magnitudes and several flares > 1 .0 magnitudes in each M-band exposure. We can also 
predict the probability of observing a large flare occurring on a star that is too faint 
to be seen during quiescence. These large flares, which appear as optical transients in 
the survey, have a few percent probability of being observed in each M-band exposure, 
although the uncertainties are large. 

Finally, when large numbers of flares are measured in time domain surveys such 
as Pan-STARRS, PTF, and LSST, our model can be used to precisely determine the 
FFD of stars in very small bins of spectral type or color and distance from the Galactic 
plane. This will inform our knowledge of how the magnetic field of M dwarfs changes 
with spectral type (in particular how it changes across the fully convective boundary), 
with age, and with activity level. 
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